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We present quantum Maxwell-Bloch equations (QMBE) for spatially inhomogeneous semicon- 
ductor laser devices. The QMBE are derived from fully quantum mechanical operator dynamics 
describing the interaction of the light field with the quantum states of the electrons and the holes 
near the band gap. By taking into account field-field correlations and field-dipole correlations, the 
QMBE include quantum noise effects which cause spontaneous emission and amplified spontaneous 
emission. In particular, the source of spontaneous emission is obtained by factorizing the dipole- 
dipole correlations into a product of electron and hole densities. The QMBE are formulated for 
general devices, for edge emitting lasers and for vertical cavity surface emitting lasers, providing a 
starting point for the detailed analysis of spatial coherence in the near field and far field patterns of 
such laser diodes. Analytical expressions are given for the spectra of gain and spontaneous emission 
described by the QMBE. These results are applied to the case of a broad area laser, for which the 
frequency and carrier density dependent spontaneous emission factor /3 and the evolution of the far 
field pattern near threshold are derived. 

PACS numbers: 42.55.Px, 42.50.Lc 



I. INTRODUCTION 

The spatio-temporal dynamics of semiconductor lasers can be simulated successfully by semiclassical Maxwell-Bloch 
equations without including any quantum effects in the light field (for an overview of the theory and modeling see 
and references therein). The classical treatment of the light field is justified by the high intensity of the laser light 
well above threshold. However, the incoherent noise required by the uncertainty principle in both the electrical dipolc 
of the semiconductor medium and the light field itself is of significant importance when several cavity modes compete 
or when the laser is close to threshold. 

Photon rate equations for multi-mode operation of semiconductor lasers show that the spontaneous emission terms 
may contribute significantly to the spectral characteristics of the light field emitted by the laser (2||j|. Such models 
assume a fixed mode structure determined entirely by the empty cavity. This assumption does not apply to gain 
guided lasers and to unstable resonators, however In these cases it is therefore desirable to explicitly describe 

the spatial coherence of spontaneous emission. 

The spatial coherence of spontaneous emission and amplified spontaneous emission is even more important in devices 
close to threshold or devices with a light field output dominated by spontaneous emission such as superluminescent 
diodes and ultra low threshold semiconductor lasers . Ultra-low threshold lasers may actually operate in a regime 
of negative gain where spontaneous emission is the only source of radiation . A description of the light field emitted 
by such devices therefore requires an explicit description of the spatial coherence in spontaneous emission as well. 

An approach to the consistent inclusion of the quantum noise properties of the light field in the dynamics of 
semiconductor laser diodes using nonequilibrium Green's functions has been presented in |p[-pj]|. In this approach, 
the linear optical response of the medium is varied as a function of the time-dependent electron-hole distributions. 
Although the non-equilibrium Green's function presents an elegant solution for the description of many-body effects 
|^| , the representation of the interband dipole dynamics by Green's functions causes a non-Markovian memory effect 
which is difficult to handle and is therefore usually neglected JlOfl , Moreover, the need to determine the Green's function 
corresponding to the dynamically varying carrier distribution requires a computational effort far greater than that 
required for the integration of the corresponding Maxwell-Bloch equations. Therefore, as stated in pJj ], an exact 
analytical investigation of the spatial mode structures in realistic cavities using non-equilibrium Green's functions is 
out of reach. In order to simulate the spatiotemporal dynamics of multi mode operation, of lasers near threshold, 
low threshold lasers or superluminescent diodes, it is therefore dcsircablc to formulate an alternative approach to the 
problem of spontaneous emission and amplified spontaneous emission in such devices which is based on Maxwell-Bloch 
equations. By including the spatiotemporal dynamics of the interband dipole in such equations, non-Markovian terms 
are avoided and the quantum mechanical equations may be integrated in a straightforward manner. 



- * — < 

X 

5-H 



1 



The starting point for our description of quantum noise effects is the dynamics of quantum mechanical operators of 
the field and carrier system. Since the operator dynamics of the carrier system have been investigated in the context 
of Maxwell-Bloch equations before |0| and the light field equations correspond exactly to the classical Maxwell's 
equations, it is possible to focus only on the local light-matter interaction. Once the properties of this interaction 
are formulated in terms of the expectation values of field-field correlations, dipole-field correlations, carrier densities, 
fields and dipoles, the dynamics of the carrier system and the light field propagation may be added. 

In section [n| the quantum dynamics of the interaction between the light field and the carrier system is formulated in 
terms of Wigner distributions for the carriers and of spatially continuous amplitudes for the light field. The equations 
are formulated for both bulk material and for quantum wells including the effects of anisotropic coupling to the 
polarization components of the light field. Section |l| summarizes the effects of the dynamics of the electron-hole 
system in the semiconductor material. The light field dynamics are introduced in section IV . By quantizing Maxwell's 
equation, the coupling constant go introduced in section || is expressed in terms of the interband dipolc matrix element. 
The complete set of quantum Maxwell-Bloch equations is presented in section [y| Based on this general formulation, 
specific approximate versions for quantum well edge emitting and vertical cavity surface emitting lasers are derived. 
The possibility of including two time correlations in the quantum Maxwell-Bloch equations is discussed and equations 
are given for the case of vertical cavity surface emitting lasers. In section |Vl| , analytical results for the spectra of gain 
and spontaneous emission in quantum wells as well as the spontaneous emission factor (3 and the far field pattern of 
amplified spontaneous emission in broad area quantum well lasers are presented. Section VII concludes the article. 



II. DYNAMICS OF THE LIGHT-CARRIER INTERACTION 



A. Hamiltonian dynamics of densities and fields 

In following, we will describe the active semiconductor medium in terms of an isotropic two-band model where, for 
the case of the holes, a suitably averaged effective mass is taken [Q . A generalization to more bands is straightforward. 
In terms of the local annihilation operators for photons (6r), electrons (cr), and holes (dn), the Hamiltonian of the 
light-carrier interaction can be written as 
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The operator dynamics associated with this Hamiltonian are then given by 
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The discrete positions R correspond to the lattice sites of the Bravais lattice describing the semiconductor crystal. 
Each lattice site actually represents the spatial volume vq of the Wigner-Seitz cell of the lattice. For zincblende crystal 
structure, this volume is equal to one quarter of the cubed lattice constant. In the case of AlGaAs structures the lattice 
constant is about 5.65 x I0~ 10 m and w 4.5 x 10~ 29 m 3 |l3). The photon annihilation operator S R therefore describes 
the annihilation of a photon within a volume v§. Here, we will focus on the light-carrier interaction and the quantum 
noise contributions responsible for spontaneous emission. For that purpose, we extend the semiclassical description 
by including not only the expectation values of the field and dipole operators, (b R ) and (c R <i R /), respectively, but also 
the field-field correlations (^r&r') and the field-dipole correlation (6 R c R 'd R »). The factorized equations of motion 
then read 
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Note that this set of equations already represents a closed description of the field dynamics. If, as in many experimental 
configurations, the absolute phase of the light field and dipole operators may be considered unknown, these equations 
are sufficient for a description of the light-carrier interaction. However, when two time correlations are of interest or 
in the case of coherent excitation by injection of an external laser it may also be necessary to additionally consider 
the dynamics of the field and dipole expectation values, i.e., 
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B. Physical background of the factorization 

In the following, we will briefly discuss the implications of the factorization performed in the derivation of (|3^). 
In order to formulate the dynamics of the light-matter interaction without including higher order correlations, three 
terms have been factorized in the time derivative of the field-dipole correlation (^). The factorizations are 

(6 R 6 R » c R „ c R < ) S3 (6 R 6 R » ) (c R » c R < ) (4a) 

(b^bR, d\, 4" ) « {btJ>Ri ) (4 d R" ) ( 4b ) 

(c R c R /d R ^R"> ~ (c R 5 R')(4^R")- ( 4c ) 

No additional factorizations are necessary in the density dynamics of photons, electrons and holes. The three fac- 
torizations are based on the assumption of statistical independence between the densities of photons, electrons and 
holes. Since we are not considering fluctuations in the particle densities, this is a necessary assumption. 

Equations ( f4a| ) and ([it]) separate the photon density from the carrier densities. These terms represent stimulated 
emission processes. Therefore, a correlation of the photon density with the carrier densities would lead to a modified 
stimulated emission rate. Below threshold, this effect will be small because the amplified spontaneous emission is 
distributed over many modes such that the local correlations between photon and carrier densities are weak. Above 
threshold, the photon number fluctuations in the lasing mode cause relaxation oscillations. The photon number 
fluctuations are nearly ninety degrees out of phase with the carrier number fluctuations. Therefore, the time averaged 
correlation is still negligible. 

Equation (^c]) separates the electron and hole densities. This term represents the spontaneous emission caused by 
the simultaneous presence of electrons and holes in the same location. Although it is reasonable to assume that the 
high rate of scattering at high carrier densities effectively reduce all electron-hole correlations to zero, it is important 
to note that the interband dipole (c R o£ R ') implies a phase correlation between the electrons and the holes. In fact, 
the spontaneous emission term factorized according to equation (^) originates from the dipole-dipole correlation 
(c R 4'^R"^R"')- Note that this term could also be factorized into the product of dipole operators (c R d R /)*(c R »d R <»). 
The dynamics of the field- field and the field-dipole correlations are then identical to the dynamics of the products of the 
fields and dipoles. Therefore, that factorization corresponds to the approximations of the conventional Maxwell-Bloch 
equations such as described in Jl2| which do not rigorously include spontaneous emission. 

Generally, spontaneous emission must always arise from random phase fluctuations. These are given by the product 
of electron and hole densities. While the phase dependent dipole relaxes quickly due to scattering events, the carrier 
densities are preserved during scattering. Therefore 

(c R d R ')*(c R «d R /») < (c R CR")( rf R 'dR"') (5a) 

is usually a good assumption in semiconductor systems. Note that this assumption does fail in the case of low carrier 
densities and high dipole inducing fields. However, this case only occurs if the light field is injected from an external 
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source. In semiconductor lasers and in light emitting diodes, the major contribution to the dipole-dipole correlations 
stems from the product of electron and hole densities. To check the statistical independence of photon, electron, and 
hole densities, it is also convenient to check the corresponding inequality for the three particle coherence represented 
by the field-dipole correlation, 

(ty l fa>dw>)*{ty l mfo>'>d R '>>>>) < (&R&R»')(4' 5 R''' , )(4" ( fe'''''>- (5b) 

Thus, if a calculation does not fulfill this requirement, particle density correlations additionally have to be taken into 
account. 



C. Wigner function formulation 

In order to connect the light-carrier interaction to the highly dissipative carrier transport equations, it is practical to 
transform the carrier and dipole densities using Wigner transformations ]l4| . Replacing the discrete density matrices 
by continuous ones obtained by polynomial interpolation will allow e.g. for numerical purposes an arbitrary choice 
of the discretization scales, which generally will be much larger than a lattice constant. Analytically, it permits an 
application of differential operators. Physically, the particle densities are smooth functions over distances of several 
lattice constants. A coherence length shorter than e.g. ten lattice constants would require k states with | k | of at least 
one twentieth of the Brilloin zone diameter. In typical laser devices, however, the electrons and holes all accumulate 
near the fundamental gap at k = 0. Therefore, it is simply a matter of convenience to define the continuous densities 
such that 

— (4cr.'} (6a) 
-<4d Rf ) (6b) 

— (cr<2r/). (6c) 
These continuous functions may then be transformed into Wigner functions by 
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f h (r, k) = J dVe^'/ ( r ~ r + ^) ( 7b ) 
p(r, k) = J d 3 r'e-^'p d ^ le (r - -, r + -V (7c) 

The normalization of these Wigner functions has been chosen in such a way that a value of one represents the 
maximal phase space density possible for Fermions, that is one particle per state. Since the density of states in the 
six-dimensional phase space given by r and k is 1/87T 3 , a factor of 1/87T 3 will appear whenever actual carrier densities 
need to be obtained from the Wigner functions. However, the normalization in terms of the maximal possible phase 
space density is convenient because it represents the probability that a quantum state in a given region of phase space 
is occupied. Therefore, the Wigner distribution corresponding to the thermal equilibrium of a given particle density 
is directly given by the Fermi function. 

To deal with the light field dynamics in the same manner, the field and field-field correlation variables must also be 
defined on a continuous length scale. In order to obtain photon densities, we define 

£(r = R) = -^(b R ) (8a) 



J(r = R;r' = R') = -<6 R 6R'>. (8b) 
Finally, the dipole-field correlation must be defined accordingly, such that 

& corr (r - R;r' - R',r" = R") = -L^cr^r") (9) 

C(r; r', k) = J d 3 r" e - ikr " Q corr - (r; r' - y , r' + y j . (10) 
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With these new definitions, the light-carrier interaction dynamics can now be expressed in a form which considers 
both the position and the momentum of the electrons and holes. The dynamics of emission and absorption now reads 
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D. Local approximation 



The integrals over x and q represent seemingly non-local effects introduced by the transformation into Wigncr 
functions. This property of the Wigner transformation retains the coherent effects in the carrier system. For the 
interaction of the carriers with the light field, it ensures momentum conservation by introducing a non-local phase 
correlation in the dipolc field corresponding to the total momentum of the electron and hole concentrations involved. 
Effectively, the integral over q converts the momentum part of the Wigner distributions into a coherence length. This 
coherence length then reappears in the spatial structure of the dipole field and the electromagnetic field generated by 
the carrier distribution. However, the coherence length in the carrier system is usually much shorter than the optical 
wavelength. It can therefore be approximated by a spatial delta function. Here, we do this by noting that 
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If the effects of the momentum shift q in the Wigner functions is neglected, the integrals may then be solved, yielding 
only local interactions between the carrier system and the light field: 
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These equations now provide a compact description of the light-carrier interaction in a three dimensional semicon- 
ductor medium, including the incoherent quantum noise term which is the source of spontaneous emission. 



E. Light-carrier interaction for quantum wells 

Similar equations may also be formulated for a quantum well structure by replacing the phase space density of 1/8-7T 3 
with 1/47T 2 , reducing the spatial coordinates of the carrier system to two dimensions, and introducing a delta function 
for the coordinate perpendicular to the quantum well at the points where field coordinates correspond to dipolc 
coordinates. Of course, the electromagnetic field remains three dimensional, even though the dipole it originates from 
is confined to two dimensions. In particular, the correlation C(r; r', k) has both a three dimensional coordinate r and 
a two dimensional coordinate r'. It is therefore useful to distinguish the two dimensional and the three dimensional 
coordinates. In the following, the two dimensional carrier coordinates will be marked with the index ||. Note that in 
some cases, both r and rn appear in the equations. In those cases, the in-plane coordinates r x and r y are equal, while 
the perpendicular coordinate r z must be equal to the quantum well coordinate Zq. The equations for the interaction 
of the three dimensional light field with the two dimensional electron-hole system in a single quantum well subband 
then read 
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Note that the value of go will usually be slightly lower than the bulk value because the overlap of the spatial wave- 
functions of the electrons and the holes in the lowest subbands is less than one. The equations derived above represent 
the interaction of a single conduction band and a single valence band with a single scalar light field. Neither the spin 
degeneracy of the carriers nor the polarization of the light field has been considered. 
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F. Spin degeneracy and light field polarization 



Since the geometry of light field emission is highly dependent on polarization effects such effects should also be 
taken into account in the framework of this theory. The basic interaction between a single conduction band, a single 



valence band and a single light field polarization are accurately represented by equations (13a -131) and (14a -141). By 
adding the contributions of separate transitions, any many band system may be described based on these equations. 
In semiconductor quantum wells the situation is considerably simplified if only the lowest subbands are considered. 
Then there are only two completely separate transitions involving circular light field polarizations coupled to a single 
one of the two electron and hole bands. The quantum well structure does not interact with light fields which are 
linearly polarized in the direction perpendicular to the plane of the quantum well. The equations for quantum wells 
are therefore completed by adding an index of + or — to each variable. 

The situation in the bulk system is much more involved. The transitions occur between the two fold degenerate 
spin 1/2 system of the electrons and the four fold degenerate spin 3/2 system of the holes. All three polarization 
directions of the light field are equally possible, connecting each of the electron bands with three of the four hole 
bands. However, since the effective mass of the two heavy hole bands is much larger than the effective mass of the 
light holes (e.g. by a factor of eight in GaAs), only a small fraction of the holes will be in the light- hole bands (about 
6% in GaAs for equilibrium distributions). Consequently, the carrier subsystem can again be separated into two 
pairs of bands. However, the light field polarization emitted by the electron-heavy hole transitions in bulk material is 
circularly polarized with respect to the relative momentum 2k of the electron and the heavy hole. Since usually there 
is no strong directional anisotropy in the k space distribution of the carriers, it can be assumed that one third of the 
k space volume contributes to each polarization direction and the equations may be formulated accordingly. 

In the following, we will assume that the Wigner distributions of the two pairs of bands considered are approximately 
equal at all times. Note that this means that hole burning effects in the spin and polarization dynamics which may 
occur in vertical cavity surface emitting lasers [fl5| are ignored. However, such effects have been investigated in several 
other studies |l6|-|l9| and are found to be fairly weak in some devices |2(J . 

The complete set of dynamical equations can now be formulated by adding the carrier dynamics and the linear part 
of Maxwell's equations to the light-carrier interaction. 

III. CARRIER DYNAMICS 

Modeling the carrier dynamics of a semiconductor system can be a formidable task all by itself. A number of 
approximations and models have been developed to deal with the effects of many-particle interactions and correlations 
and with the dissipation caused by the electron- phonon interactions |2l],|2^]. In the following, we choose a simple 
diffusion model. Many particle effects such as the band gap renormalization or the Coulomb enhancement are not 
mentioned explicitly, but can be added in a straightforward manner fUM. 

We assume that the electron and hole densities will be kept equal by the Coulomb interaction, which will induce a 
current whenever charges are separated. Therefore, it is possible to define the ambipolar carrier density -/V(r) with 

^) = i/rf 3 k/ e (r,k) 



^LJd 3 kf h (r,k). (15) 



Note that the two fold degeneracy of the electron and heavy hole bands has been included by choosing a density of 
states of j^j instead of -^j. This includes the assumption that the Wigner distribution does not depend on the spin 
variable of the electrons and holes as mentioned above. 



The light carrier interaction of this carrier density is 
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where the index i denotes the component of the light field or dipole density corresponding to the linear polarization 
direction of i = x,y,z. In the case of CV/(r;r',k) the index i refers to the field polarization and the second index 
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j denotes the vector component of the dipole vector. Equation (|l^) shows how the field-dipole correlation converts 
electron- hole pairs into photons. The total carrier density dynamics can now be formulated as 



D amb AN(r)+j(r)-jN(r 



47T 3 



J d 3 k^(C»(r;r,k)-C^(r;r,k)), 



(17) 



where D am i> is the ambipolar diffusion constant, j(r) is the injection current density, and 7 is the rate of spontaneous 
recombinations by non-radiative processes and/or spontaneous emission into modes not considered in iy(r, r'), e.g. 
if the paraxial approximation is applied. 

The k dependence of the distribution functions / e (r, k) and /' l (r, k) may be approximated by assuming that the 
electrons and holes will always be in thermal equilibrium. The distribution functions are then given by Fermi functions 



/£ fc (r,k)= exp 



1 



k B T \ 2m 



e . ll 
eff 




(18) 



where rn^j are the effective masses of electrons and heavy holes, respectively. The chemical potential fi e ' (r) is a 
function of the carrier density N(r). A useful estimate of this relationship is given by the Pade approximation |p"2| , p3f . 
Spectral holeburning may be taken into account by introducing a relaxation time r r and converting the dynamics 
of the distribution function due to the light-carrier interaction into a deviation from the equilibrium distribution by 
adiabatic elimination of the relaxation dynamics: 



/^(r,k)=/ e ^(r,k) 



.£(a(r;r,±k)-C*(r;r,±k)). 



(19) 



Finally, the carrier dynamics of the dipole p(r, k) and the dipole part of the field-dipole correlation C(r; r', k) needs 
to be formulated. Since both depend on a correlation of the electrons with the holes, they will necessarily relax 
rather quickly at a rate of T(k) which should be of the same order of magnitude as l/r r . Physically, T(k) may be 
interpreted as the total momentum dependent scattering rate in the carrier system. The remainder of the dynamics 
can be derived from the single particle dynamics. This unitary contribution to the evolution of the dipole may be 

expressed by a momentum dependent frequency Q, (k). For parabolic bands and isotropic effective masses rrfjh, this 
frequency is given by 
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Many-particle effects due to the Coulomb interaction between the carriers may be included by introducing a carrier 
density dependence in r(k, 7V(r)) and f2(k, N(r)). Such renormalization terms representing the mean field effects 
of the carrier-carrier interaction have been derived and discussed e.g. in [jl2j. In the following, this many particle 
renormalization will not be mentioned explicitly, although it can be included in a straightforward manner. 
With the rates T(k) and f2(k) the dipole dynamics reads 



d_ 

at 



C 4J (r;r',k) 
J^(r,k) 



= -(r(k) 



= -(T(k) 



in(k))Cy(r;r',k) 
ifi(k))pi(r,k). 



(21a) 



(21b) 



Note that the phase dynamics is formulated relative to the band gap frequency ujq. The real physical phase oscillations 
of p(r, k) would include an additional phase factor of exp[— iu>ot\. However, the only physical effect of this oscillation 
is to establish resonance with the corresponding frequency range in the electromagnetic field, the dynamics of which 
we consider next. 



IV. MAXWELL'S EQUATION 

The Heisenberg equations of motion describing the operator dynamics of the electromagnetic field operators are 
identical to the classical Maxwell's equations. In terms of the electromagnetic field E(r) and the dipole densities P(r), 
the equation reads 
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Vx(VxE(r)) + ^(E(r) + ^-P(r))=0, (22) 

where eo and c are the dielectric constant and the speed of light in vacuum, respectively, and e r eo is the dielectric 
constant in the background semiconductor medium. 

Maxwell's equation describes the light field dynamics for all frequencies. Since we are only interested in fre- 
quencies near the band gap frequency lu , it is useful to separate the phase factor of exp[— Wot], defining E(r) = 
exp[— ia;o£]Eo(r). Now Eo(r) can be considered to vary slowly in time relative to exp[— iuot]. Therefore, the time 
derivatives may be approximated by 

d 2 d 
exp[iiv t]— exp[-iuj i\E (r) » -u> E (r) - 2iw — E (r). (23) 

Similarly, Po(r) may be defined such that P(r) = cxp[— itJo^]Po( r )- The approximation used here may even be of 
zero order, since we are primarily interested in the dynamics of the electromagnetic field: 

d 2 

exp[iuj t}-^ exp[-iuj t]P {r) « -uj 2 P {r). (24) 
The temporal evolution of the electromagnetic field now reads 

^Eo(r) = -i^- (V x (V x E (r)) - e r fc 2 E (r)) - i-^P (r), (25) 



where fco = ujq/c is the vacuum wavevector length corresponding to ujq. In (|25|), the field dynamics is described in 
terms of electromagnetic units, that is the fields represent forces acting on charges. To switch scales to the photon 
densities represented by £(r), energy densities have to be considered. Since the energy of each photon will be close 
to the bandgap energy fkoo, the energy density of the electromagnetic field is given by 

M)f*(r)f(r) = ^ES(r)E5(r). (26) 

Therefore, the field may be expressed as photon density amplitude using 



Eo(r) = J 2 -^£(r). (27) 

V £r£0 

The dipole density Po( r ) may be expressed in terms of p dl P°' e (r, r) and p(r, k) by noting that the density p dl P° le (j^ r ) 
is the dipole density in units of one- half the atomic dipole given by the interband dipole matrix element d cu at k = 0. 
The factor of one-half is a logical consequence of the property that (cr^r/) < 1/2. Thus, a fully polarized lattice 
would have a dipole density of / r, dl P° Je (r, r = 1)/(2vq) which must correspond to P(r) = & cv /v$. Note that d cv contains 
an arbitrary phase factor depending only on the definition of the states used for its determination. For convenience, 
we assume a definition of phases such that d cv is real. The dipole density Po(r) may then be written as 



P(r) = 2d cv p dlpole 
2 | d cv 



8tt 3 

Written in terms of £ (r) and p(r, k), the complete field dynamics now reads 



J d 3 kp(r,k). (28) 



= -i^T (V x (V x S(r)) - irkt£(T)) - i 9o |5 J d'k p(r,k). (29) 



The coupling frequency go introduced in equation(fij) may be expressed in terms of the dipole matrix element d c 



.90 = «/ "° | dcv | ■ (30) 

With this equation, the operator dynamics of the light field operator 6r corresponds to the field dynamics of the 
Maxwell-Bloch equations for classical fields. By applying the linear propagation dynamics of the field to %(r; r') and 
Cij(r; r', k) as well, it is now possible to formulate a complete set of quantum Maxwell-Bloch equations. 
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V. QUANTUM MAXWELL-BLOCH EQUATIONS 



A. Quantum Maxwell-Bloch equations for a three dimensional gain medium 

On the basis of the discussion in the previous sections, the quantum Maxwell-Bloch equations for a bulk gain 
medium in three dimensions read 

9 N(t) = D amb AN(r) + j(r) - 7 7V(r) 



Ot 



^8 



J d 3 k£(C^r;r,k)-CS(r;r,k)) (31a) 



JWr; r', k) = - (r(k) + iO(k)) C (r; r', k) 



at 



E a^7 £r 1 a^7 c ' l ^ r;r '' k ) 



S w^Wh ^ r '' k) + fc ° 2Cu (r; r '' k) 



+*2o^ (f e eq (A; iV(r')) + /* (fc; iV(r')) - l) (r; r') 



+*.9o^ <5(r - r')%/ e e 9 (k; N(r)) ■ /* (k; N(r)) (31b) 



g | d 8 k(Cy(r;r',k)-C; i (r';r,k)) (31c) 



-.9o- 87 



^p i (r > k)=-(r(k)+ifi(k))p i (r > k) 



(fl g (k; N(r)) + /* (fc; 7V(r)) - l) £(r) (31d) 

-^o|^ / d 3 k Pi (r,k). (31e) 

In order to describe a realistic diode one needs to describe not only the volume of the active region but also the 
propagation of light out of this region. This may either be achieved by defining realistic boundary conditions or by 
setting all material properties equal to zero outside a finite active volume and calculating the light field propagation 
into the outside medium by varying e r in space. 

B. Three dimensional quantum Maxwell-Bloch equations for quantum wells 

Next, we will formulate the equations for a quantum well structure. For this case, we will also describe different 
cavity structures and the appr opriate paraxial approximations possible for the various types of laser devices. Using 



the terminology of equations ( 14a -141), the quantum Maxwell-Bloch equations for a multi quantum well structure 
with Q quantum wells read 
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^iV(r||) = D amfe AiV(r||) +i(r,|) - 7 iV(r|| 



/j/ 



J d2k \\ E ( c «ii( r ; r ii ^11)^=20 - c *i||( r ; r ih k nk=*o) ( 32a ) 



9 C y|| (r;r'||,k||) = - (r(k„) + ifi(k|,)) C y -„ (r; r'y , k„) 



dt 



- E ^^C kH (r; r'|| ,k„) + fcgC^ (r; r' „ , k„)^ 



dt ij[r ' r) ~ l 2kl^\dr k tr 8r k dr' k er dr'J 



d 

( r ll > k ll ) = " ( r ( k ll ) + ^( k ll )) Pi\\ ( r ll ' k ll ) 



dt 



(32b) 



■ ^0 V ( 9 -1 9 r / / \ 9 -I 9 t ( l\\ 

' l9Qy § J d2k " (^ii( r ; r V k ii)<^-zo) 

-C; il| (r';r||,k,|)5(r,-«b)) (32c) 



/" d 2 k|, Pi|| (r||,k||)5(r z -0 O )- (32e) 



Note that iV(r||) is the total carrier density. Therefore, the density per quantum well which determines the chemical 
potential of the carrier distribution functions is iV(r|| )/Q. Again, the structure of an external cavity may be considered 
either by boundary conditions or by spatially varying e r . In particular laser diodes may be described by distinguishing 
between reflective and non-reflective edges. If the reflective surface is perpendicular to the plane of the quantum wells, 
the laser is an edge emitter. If the reflectivity is very high on the surface planes parallel to the quantum wells, the 
laser is a vertical cavity surface emitting laser (VCSEL). 



C. One dimensional quantum Maxwell-Bloch equations for edge emitting lasers 

In an edge emitting laser, the laser light field propagates in the plane of the quantum well. Since the z-direction 
is already defined as the one perpendicular to the quantum well we will define the y-direction as the direction along 
which the laser light propagates. A schematic representation of this type of laser geometry is shown in figure [I]. It is 
possible to drastically reduce the dimensionality of the equation describing the edge emitting laser geometry by noting 
that the light field polarization of the amplified fields will be in the plane of the quantum well and by limiting the 
analysis to a single longitudinal mode. Effectively, this corresponds to a light field £(r) with the following properties: 
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£ x (r) := £o(r x )£(r y , r z ) (33a) 
£ z (r):=0 (33b) 

The envelope function £(r y ,r z ) describes both the propagation along the y-direction and the confinement along the 
z-direction. It represents an approximate solution of the wave equation in the yz-plane normalized by 



/ 



dr y dr z | £(r y ,r z ) | 2 = 1. (33d) 



The equations are then limited to light field modes with the two dimensional envelope £(r y , r z ). Spontaneous emission 
into other light field modes must be considered by including the rate of emission in the carrier recombination rate 
7. Since the length L of the laser in the y-direction is also an important property of the device, it is included by 
considering the openness of the optical cavity. With the reflectivities of the laser mirrors given by R\ and R2, the 
light field in the cavity is damped by losses through the mirrors at a rate of 

K = -^=ln[fliifc]. (34) 
The new one-dimensional variables are now defined as follows: 

(35a) 
(35b) 
(35c) 
(35d) 
(35e) 

The carrier density is now given in terms of a one dimensional density. To obtain the two dimensional carrier density 
per quantum well, this density is to be divided by QL. Note that the intensity is also given in terms of photons per 
unit length. The dynamics of the edge emitter then reads 

d d 2 

-^N 1D (r x ) = D amb -^N 1D (r x ) + Lj(r x ) - -fN 1D (r x ) 

+ig °^~ / d2k|1 ( C o( r x\r x ,^\\) ~ c o( r x\r x ,iz\\)) (36a) 
^C (r x ;r;,k||) = - (r(k„) + ifi(k|,)) C (r x ; r' x , k„) 



N 1D (r x ) = 


1 dryN^) 




C (r x ;r£.,k||) = ^ 


[ dr y dr z dr' v £(r y ,r z )C{r v ,r' z 


= 2;o)C xx (r;r'||,k||) 


h{r x ;r' x ) = 


[ dr y dr z dr' y dr' z £(r y ,r z )C(r y 


,r' z )I xx (r;r') 


Po(r-x,k||) = ^ 


[ dr y C(r v ,r z = z )Px(r||,k||) 




£o(r x ) = 


/ dr y dr z C(r y ,r z )£ x (r). 





dt 



up d 2 
^k 2 e r dr 2 x 



-/gCb(rx;ri,k||) - i —j- ^ C (r x ; r' x , k,,) 



(A (*„; ^) + /* (*„; ^) - l) Io(r,;0 



Nip(r x ) \ h ( N lD {r x ) 
QL )-t«i\ k h ql 



+ig a^T S(r x - r' x )f eq [k\\ ; • ( fc„ ; ) (36b) 



Wt I»{r x y x ) = - 2K I (r x y x )-i^- - _j / o(r ,y x) 

-fco^j I d 2 k|, (C (r x ;r;,k||)-C *(r x ;r x ,k||)) (36c) 



5 

— p (r x ,k||) = - (r(k||) + in(k||))po(rx,k|| 
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+<W«* ( /«, *||5 + /* fc||; ) - 1 ) £b(r.) (36d) 



Nip(r x ) \ h ( N lD {r x ) 

^. £ o{r x ) = -KE {r x )+i-^^^Eo{r x )-igQ^^ ; J d 2 k|| po{r x ,'k\\), (36e) 
where a is the confinement factor which determines the overlap between the quantum wells and the light field mode, 

<j = Q J dr y \ £(r y ,r z = z ) | 2 . (37) 



D. Two dimensional quantum Maxwell-Bloch equations for VCSELs 

In a VCSEL the light field is strongly confined by highly reflective mirrors above and below the quantum wells. The 
light field propagates perpendicular to the quantum well structure as shown in figure ||. Therefore both the possible 
polarization directions and the spatial dynamics remain two dimensional. Only the ^-direction may be eliminated by 
averaging over a single longitudinal mode. Coupling terms between the polarization directions should be taken into 
account, even if they are small. For VCSELs, the assumptions read 

£{v) « £(r u )£(r,) (38a) 

le z (r) « - A [e x{r) + ^ / d%si r, - , K (r||)k||) ) 

-A ^(r) + J d^S(r z - Zb)p v (r||,k„)) . (38b) 

The latter condition takes into account the divergence of £(r) caused by the polarization p||(r||,k||). This is an 
important contribution to the quantum Maxwell-Bloch equations, since it coherently couples orthogonal polarizations. 
The properties of the envelope function £(r z ) are defined as 

dr z | £(»•») | 2 = 1 (38c) 

« ~k 2 ar z ). (38d) 

The cavity loss rate k is defined as in equation (Q). However, the experimentally observed polarization stability 
is taken into account by using slightly different reflectivities for the x and y polarizations. The cavity loss rate is 
therefore given by n x and K y . Experimental results [|o| suggest that (k v — n x )/ (k x ) w — 10 -2 . A birefringence of 
Su) x / y is also included to denote the difference between the band gap frequency and the longitudinal frequencies of the 
confined light field for the two polarization directions. The birefringence 5oj x — 8uj y is usually in the GHz range. Since 
all coordinates are two dimensional, the index || which marked the two dimensional coordinates in the quantum well 
equations will be omitted. Instead, the two dimensional variables are marked with a tilde. All indices and coordinates 
of such variables are only defined in two dimensions. The variables are defined as 

N{t)=N{t\\) (39a) 

C^-(r;r',k) = j dr z ^r z )C(r' z = z )C ih (r; r'y, k„) (39b) 

4-(r; r') = J dr z dr^(r z )C (r z )I^(r; r') (39c) 
Pi(r,k) =£*(r z = Zo)Pi||(r||,k||) (39d) 
£i(r) = f dr z C(r z )£ t (r). (39e) 



Note that the dipole variable Pi„ (ry, k||) is rescaled by the field density of the envelope function £(r z ) at r z — zq. The 
dipole given by f>i{v, k) is therefore the average dipole density within the whole cavity, not just within the quantum 
wells. The quantum Maxwell-Bloch equations for the VCSEL now read 
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^iV(r) = D amb AN(r) +j(r) - jN(r) 



9 (r; r', k) = - (T(k) + ifi(k)) <5y(r; r', k) 



- Wwj)Cij (r;r',k) - i— 5— A r CV,-(r;r',k) 



+i9ovV^ [fl q U ^p-j + \k; ^f-j - lj 4(r; r') 

S(t - r')%/ e e g [A; ^ J • ( A; ^ J (40b) 



Q 7 T Q 



— Jjj ;(r; r') = + Kj)!^; r' - (A r - A r /) 7 y (r; r') 



-^o^5 | d 2 k(Cy(r;r',k)-C; i (r';r,k) 



C lfc (r;r',k)-^— C; fc (r';r,k) (40c) 



Pi(r,k) = -(r(k)+in(k))pi(r,k) 



+iW** ^/ e e 9 \k; ^ j + / e " g (a; ^ j - lj £(r) (40d) 
jrA{*) = + iSLJi)£i(T) + i Jf° Agi(r) 



with (j being the confinement factor along the z-direction, 

0- = Q I £(r, = 20) I 2 • (41) 

Equations (|4^) present a starting point for the study of spatial polarization patterns and fluctuations in VCSELs. For 
more realistic models, it may also be desirable to include a spatial dependence of the birefringence and the dichroism. 
Also, nonlinear effects may be introduced, e.g. by separating the carrier densities for right and left circular polarization 



E. Statistical interpretation and two time correlations 



Using the equations presented above, it is now possible to calculate the emergence of a spatially coherent light field 
in semiconductor laser diodes both above and below threshold. Note that the average light field E(r) will remain zero 
at all times due to the random phases of spontaneous emission processes. The (average) spatial coherence of the light 
field, however, does not vanish and is fully described by the non-local field-field correlations I(r; r') which emerge due 
to the propagation and/or amplification of the originally incoherent local spontaneous emissions. This emergence of 
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coherence as a concequence of incoherent emissions has been discussed in a temporal context using nonequilibrium 
Green's functions in (Hj. 

In order to understand the physical implications of the well known absence of an average light field in lasers, it 
should be recalled that all the results of the quantum Maxwell-Bloch equations represent averages which have to be 
interpreted in terms of statistical physics. For example, the field-field correlation I(r; r') represents a variance of the 
probability distribution with respect to the possible spatial electromagnetic field values. The coherent field observed 
in experimental time-resolved measurements will vary randomly from measurement to measurement according to this 
probability distribution. Indeed, the intensity distribution itself will vary depending on the random phase interference 
of the eigenmodes given by I(r;r'). The calculated average spatial intensity distribution I(r;r) only describes the 
average near field pattern, which is likely to be close to but not identical with the one actually observed. The 
fluctuations of the actual intensity distribution around this average, however, are disregarded as a consequence of the 
factorization performed in section |l]. Moreover, the fluctuations in the carrier density distribution induced by spatial 
holeburning associated with these fluctuations of the intensity distribution have also been disregarded. 

While the average spatial coherence of the light field is fully described by the quantum Maxwell-Bloch equations, 
the average expected temporal coherence has not yet been explicitly considered. Indeed, it is not necessary to consider 
temporal coherence in the closed sets of quantum Maxwell-Bloch equations given above because all the information 
required to obtain the correct emission and absorption rates are incorporated in the field-dipole correlation C(r; r', k). 



As shown in section VI, the spectra of gain and spontaneous emission are implicitly given by the coherent dipole 
dynamics which enters into the temporal evolution of this field-dipole correlation. If explicit information about the 
two time correlations is desired, however, such correlations may be included in the dynamics by noting that the 
temporal evolution of the two time correlations of the field I(r, t;r',t') and the two time correlations of the field- 
dipole correlation C(r, t; r', k, t') is equivalent to the dynamics of the field and dipole expectation values plfl . While 
the quantum Maxwell-Bloch equations for the correlations at t = t' remain unchanged, the evolution of the two time 
correlations as a function of t' > t is then given by an additional pair of equations which depend on the carrier 
dynamics given by the solution of the original system of quantum Maxwell-Bloch equations as presented above. In 
the case of VCSELs these additional equations supplementing equations ([h)|) read 

^A(r, t; r', k, t') = - (r(k) + ifi(k)) C VJ (r, t; r', k, t') 



4~ '", h >"* \ K ( /;„ ( k;^l\ +f: q U^l) I | /,,(.-./:,•'./') (42a.) 



- % 9^ J d 2 kC t3 (r,t;r',k,t') 
J k K o e r or-i k 



If the carrier density changes slowly, the equations describe the linear response of the medium caused by the initial 
intensity distribution t; r', t) and the initial field-dipole correlation Cy(r, t; r', k, t) determined from equations 
( fio| ) . In this case it is possible to derive the spectrum and the gain from the eigenmodes and the associated eigenvalues 
of the quasi-stationary linear optical system. Note that the eigenmodes of the linearized dynamics are not necessarily 
identical with the eigenmodes of the intensity distribution Jy (r, t;r',i), since fast variations in the carrier density 
distribution may have induced phase locking between the dynamical eigenmodes. 

In general, the frequency spectra of the light field are given by Fourier transforms of the two-time correlations 
obtained from equations (|4^). However, as noted above, such spectra do not comprise any effects which arise from 
carrier density fluctuations. These effects are known to be quite significant. Well known examples of carrier fluctu- 
ation effects in the frequency spectrum of semiconductor lasers are the linewidth enhancement phcnomcnologically 
described by the linewidth enhancement factor a and the relaxation oscillation sidebands observed in stable single 
mode operation. Moreover, spatial carrier density fluctuations may also significantly modify the multi mode spectra 
of semiconductor laser devices, as pointed out recently for the case of semiconductor laser arrays pi|. 



VI. AMPLIFIED SPONTANEOUS EMISSION PROPERTIES: ANALYTICAL RESULTS OF THE 

QUANTUM MAXWELL-BLOCH EQUATIONS 
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A. Gain and spontaneous emission 



For a given carrier distribution f e ^ h (k), the dynamics of the optical field £ (r) and the dipole density p(r) are linear. 
In this case, it is possible to integrate the equations of motion to obtain a Green's function for the field dynamics. 
The equation for bulk material reads 



= gl^- /d 3 k / dr e -(r(k)+ifi(k))r 

cL 127r 



x(f e (k) + f h (k)-l)£(r,t-T). (43a) 
Correspondingly, the equation for a multi quantum well structure of Q quantum wells has the form 

9 £{v,t) = g l9^8(r z -z ) [ d 2 k„ /°° dr e -(rCk„)-Hn(k„))r 



dt 



x (/ e (fc||) + /' l (fc||)-l)^(r,t-r). (43b) 



An expression for the rate G(u>) at which a light field mode of frequency u> is amplified can be derived by solving 
the integral over r using £ (r, t — r) « e luJT £(r, t). The real part of the result is the gain spectrum given in terms of 
amplification per unit time, G(u>). For bulk material, this amplification rate is given by 

G „ lM . ppJ^ (/»«=) + /*<*> - 1) . (44.) 

and for quantum wells, the corresponding amplification rate reads 

GqwH = 9 2 ^S(r z - z ) J d 2 k„ r2(k||) _ ^ ^h) + f h (h) - l) . (44b) 

The gain per unit length can be obtained by dividing the rate G(u>) by the speed of light in the semiconductor 

— 1/2 

medium, ce r . However, in order to establish the connection between the gain spectrum and the spectral density 
of spontaneous emission, it is more convenient to use the amplification rate as a starting point. 

The quantum Maxwcll-Bloch equations for the field-dipole correlation Cij(r;r',k) show that the ratio between the 
spontaneous contributions and the stimulated contributions is 

m C v( r '> r '> k ) I spontaneous = S(r - r%j / e (fc) ■ f h (k) 

In this equation, 5(r — r')5ij corresponds to a photon density of one photon per mode. The spectral density o f the 
spon tane ous emission may therefore be obtained by replacing (f e (k) + f h (k) — 1) with f e (k) ■ f h (k) in equations (44a) 
and ( 44t ) , respectively, and multiplying the resulting rates with twice the density of light field modes which couple to 
the medium (the factor of two being a result of considering intensities instead of fields) . At the band edge frequency 
u)q, the density of light field modes per volume and frequency interval in a continuous medium is 

uj 2 

Plight = ^Vr /2 - (46) 

The density of the spontaneous emission rate for bulk material Sbuik{oj) thus reads 

= Pu gh t 9 l^ J d 3 k r2(k) + ^ k) ^ )2 f(fc) • f h (k). (47a) 

For quantum wells, the density of modes is only 2/3 of pu g hti since the dipole component perpendicular to the 
quantum well is zero. Also, the delta function S(r z — zq) may be omitted to obtain the emission density per area. The 
spontaneous emission density Sqw{u) is then given by 

S QW {u) = 2 - PUghtg ^S(r z - zo) J d 2 k„ I ^ (k )+ r ( n^ [|) _ (<;)a / g (fe||) ' f h (h)- ( 4 7b) 
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The total rate of spontaneous emission per unit volume or area may be obtained by integrating over all frequencies. 
This integral removes the dependence on r(k). For bulk material, 

J duj S Mk (u) = pu gh t 9o-^ J d 3 k f e (k) ■ f h (k) (48a) 



and for quantum wells, 



duo S qw (cj) = PUgh t 9l^-5{r z - z Q ) \ d 2 k|, r(k [{ ) ■ f h {k^). (48b) 



For zero temperature, the spontaneous emission rate may be derived by noting that f e (k) ■ f h (k) — f e (k) — f h (k). 
The integral over k may thus be solved, resulting in 



d 3 k f e (k) ■ f h (k) = 4vr 3 iV (49a) 



for bulk material and 



J d 2 k ll t(k l] )-f h (k ]l ) = ^N (49b) 

for quantum wells. For both bulk and quantum wells the density of spontaneous emission now may be expressed as 

N 

Stotal = —, (50) 
T s 

where 1 jr s is the rate of spontaneous emission given by 

1 2?r 2 

— = -7rPUghtg v . (51) 
r s A 

Using equation (|3^), the rate of spontaneous emission may also be expressed in terms of the dipole matrix element 

dcv i 

I = J_ e V2_Li4 I dcy |2 
T s flUJo T 47T60 3c 3 CV 

4 

= T Pradi (52) 

where P ra d is the classical power radiated by an oscillating dipole of the amplitude d cv . Note that the factor of 
4/Hujo is consistent with a quantum noise interpretation of spontaneous emission such as the one represented by 
the semiclassical Langevin equations. According to this interpretation, the fluctuations of each dipole are given by 
2 | d cv | 2 because both the real and the imaginary part of the dipole contribute. The spontaneous emission of an 
excited atom is then composed of one half amplified field noise and one half dipole fluctuations. Therefore, each one of 
the two dipole components contributes one quarter of the total spontaneous emission of an excited state. Numerical 
values of | d cv | and go for GaAs may be determined by assuming a spontaneous lifetime of t s = 3 ns, a band gap of 
hu!o = 1.5 eV and e r — 12. The dipole matrix element is then | d cv |= 4.3 x 10 _29 Cto which corresponds to a distance 
of 2.7 x 10~ m times the electron charge and the coupling frequency is go = 2.1 x 10 15 s _1 . 

For quantum wells at zero temperature, the integrals over k|| can be solved analytically using equation (j2TJ) and 
assuming that T is independent of k|| . The resulting gain spectrum is then given by 



where O j is the transition frequency at the Fermi surface of the electrons and holes. It is related to the carrier density 
Nby 

ir% miff + mhff 

n / = 7? 11 S N - ^ 

1^ TO e// TO e// 
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The spectral density of spontaneous emission is thus given by 



/ s 7 2Qu m tff m eff ( /ttf-Uj\ /^\\ 

s Q wM - wo 3 ^ m e f / + 4 f/ {-^f- ) + arctan (rJ J • (55) 

Typical spectra of gain and spontaneous emission of an active layer containing Q = 5 quantum wells obtained from 
the analytic approximations (M) and (|55|), respectively, are presented in Fig. [| for characteristic values of the carrier 
density. For the spectra in Fig. H we have assumed a total spontaneous emission lifetime of r s = 3 ns and a band gap 
of %uiq = 1.5 eV. Other parameters are the effective mass of electrons (miff = 0.067mo) and holes (m^ff = 0.053mo), 
given in units of the electron mass mo as well as the dipole damping rate HT = 8 meV, and dielectric constant 
e r = 12. Gain and spontaneous emission are both displayed relative to the peak values S max — 2.7 x 10 7 cm~ 2 and 
G max — 2.5 x 10 7 cms _1 <5(r2 — zo). The gain value may be interpreted by calculating the gain of a light beam of 
width ct _1 = 10~ 5 cm with an incidence perpendicular with respect to the quantum well and traveling at a speed of 
10 10 cms~ 1 in the plane of the quantum well. The maximal gain is then given by 250cm~ 1 . The five spectra displayed 
in Fig. U are the gain at N — (go), the gain at N = 10 12 cm~ 2 (gi), the spontaneous emission at N — 10 12 cm -2 
(si), the gain at N = 5 x 10 12 cm -2 (g$) and the spontaneous emission at N = 5 x 10 12 cm -2 (s§). Fig. || clearly 
shows the influence of the carrier density N. In the absence of charge carriers (go), the laser is purely absorptive. 
With increasing carrier density N, transparency is reached at N = 10 12 cm -2 (51). This density is characterized by 
vanishing gain (i.e. there is neither gain nor absorption) at a frequency of huj = 0. At the same time, however, there 
is a significant contribution of spontaneous emission (s\) with a maximum at a frequency of hio w 4 meV. Finally at 
high values of the carrier density (N — 5 x 10 12 cm -2 ), both gain (55) and spontaneous emission (s$) have a maximum 
above the band-gap frequency. 



B. Spontaneous emission factor and far field pattern of an edge emitting laser 

In the optical cavity of a laser the equations for gain and spontaneous emission are modified by the mode structure. 
In particular, the total linear response of an electromagnetic field mode inside the cavity includes the cavity loss rate 
k. For the edge emitting semiconductor laser, the gain function of the cavity modes is given by 

GljM = s ^jA, m - ^ k[) - (A,,) + /»(,,) - 1) , (56) 

where the confinement factor a is defined according to equation (|37|). Spontaneous emission into the cavity modes 
passes through the gain medium and is thereby absorbed or amplified accordingly. Thus for an edge emitting laser, 
the rate of spontaneous emission into a cavity mode of frequency w is given by 

SM*) = tfo£*J*iH (r(k||) + ') (2k f^)^) 2 /e(fc||) - //1(fc||) - (5?) 

This rate represents the total rate of spontaneous emission events per mode, regardless of the actual width of the 
laser. On the other hand, equation (jH^) gives the total rate of spontaneous emission per quantum well area. 



1. Spontaneous emission factor 

In a laser of length L and total width W (c.f. Fig. |Tj), the spontaneous emission rate into free space is LW Stotai- 
The spontaneous emission factor (3, which is generally defined as the fraction of spontaneous emission being emitted 
into the cavity mode Q , is on the basis of our theory given by the expression 

Note that the two dimensional carrier density N in the quantum well is related to the one dimensional carrier 
density N\d by N\d = L N. The spontaneous emission factor f3 is a function of both frequency and carrier density. 
Consequently, the common assumption of the spontaneous emission factor /3 being independent of the carrier density 
may be regarded as an approximation similar to the assumption of linear gain. 
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An analytical expression for the spontaneous emission factor may be obtained for zero temperature. The k-space 
integrals may then be solved analytically using equation (|2^) and assuming that T is independent of k|| . For the 
spontaneous emission factor /?, the analytical result reads 

P(oj, Qf) = - ff TI , - _ | arctan [ -| ) + arctan [ ) ] . (59) 

,y 11 2-KpughtQ WLQ f \ \T + kJ \T + kJJ v ; 

The Fermi frequency f2/ is defined by equation (p^) and expresses, in particular, the carrier density dependence of (3. 
For f2/,(x>^;r + ftwc recover the result typically given in the literature (e.g. which does not depend on N. Fig. || 
shows the deviation of the spontaneous emission factor from this value as the Fermi frequency / passes the point of 
resonance with the cavity mode. Fig. ||| illustrates, in particular, the carrier density dependence of (3 for three modes 
with frequencies above the band gap frequency given by (a) lj = 0, (b) ui = 0.5(T + k), and (c) u> = T + k. Most 
notably, (3 is always smaller than the usual estimate given by f3(to = Q/ = 0), which is based on the assumption of 
ideal resonance between the transition frequency and the cavity mode. 



2. Far-field pattern of a broad area laser 
With equations ( |56| ) and (^), it is possible to find the steady state intensity I s of a mode with frequency oj, 

Is ^> = — — 77a- ( 60 ) 

Note that this result may also be obtained directly from equations ( |36a -36c). In a wide cavity, the cavity modes are 
approximately plane wave modes and the relation between Iq(lo) and Io(r x ; r' x ) in that case reads 

I (r x y x ) = Jdq e«*--V I (u = |^-) . (61) 

The steady state intensity distribution is characterized by the spatial coherence derived from the intensity distribution 
of the plane wave modes of the cavity. Generally, the intensity distribution of plane waves corresponds in the far field 
to an optical field at angles relative to the axis of emission in the plane of the quantum well. The angular distribution 
of intensity and coherence in the far field is thus given by 

If = / dr x dr' x e lk ^I {r Xl r' x )e-^ e '<. (62) 

Therefore, the far field intensity distribution may be determined directly from the frequency dependence of the 
intensities by 

I f (Q ) Q) = Wk I s L = ^y (63) 

the intensity is given in units of 2kTiuj$ per unit angle. 

For T — we may in analogy to ([39|) solve the integral in ([37]) by assuming V to be independent of k. As an 
analytical expression we then obtain for the far-field intensity distribution of a broad area semiconductor laser 



arctan I n ' 1 + arctan / 



1/(6,9) =Wk - 



r+K l 1 \ r+K 



tt(R + I) - 2 arctan ( %"^ B) ) 

with R = 



— arctan 



£(6) 

r+K 



', ./'./' " l eff 



glv a m e eff m^ ff 

and w(0) - (64) 

The parameter R represents the ratio between the cavity loss rate k and the maximum amplification rate of the gain 
medium. Laser activity is only possible if R < 1. The classical laser threshold is defined by the carrier density for 
which the denominator of 7/(6,0) is zero for a single specific frequency lu(Q). Consequently, the carrier density at 
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which this occurs is pinned. Figure |5| shows the far field intensity distribution for different carrier densities below this 
pinning density. In Fig. ^| (a), the wide intensity distribution of amplified spontaneous emission for carrier densities 
is much lower than the pinning density. The intensity maximum is clearly located at = 0. Figure ||(a) shows the 
intensity distribution for carrier densities halfway towards threshold. Already, the intensity maxima move to angles 
of ±15°, corresponding to the frequency at which the gain spectrum has its maximum. In the case of Fig. ^|(c), the 
threshold region is very close to the pinning density. The peaks in the far field pattern narrow as the laser intensity 
is increased. Consequently the far field pattern indeed is a measure of the spatial coherence - similar as the linewidth 
of the laser spectrum is a measure of temporal coherence. It is therefore desirable to consider quantum noise effects 
in the spatial patterns of optical systems. In the context of squeezing, such patterns have been investigated by Gatti 
and coworkers ||(| based on the general formulation of Lugiato and Castelli . The laser patterns presented here 
are based on the same principles. Usually, however, the strong dissipation prevents squeezing in laser systems unless 
the pump-noise fluctuations are suppressed pqj. 



VII. CONCLUSIONS 



The quantum Maxwell-Bloch equations (QMBE) for spatially inhomogeneous semiconductor lasers derived in this 
paper take into account the quantum mechanical nature of the light field as well as that of the carrier system. The 
only approximation used in the derivation of the intensity and correlation dynamics is that of statistical indepen- 
dence between the two carrier systems and the light field. In the QMBE presented here, the effects of coherent 
spatiotemporal quantum fluctuations which are generally not considered in the semiclassical Maxwell-Bloch equations 
for semiconductor laser devices have thus been taken into account. 

The spontaneous emission term appears side by side with the gain and absorption term in the dynamics of the field- 
dipole correlation. In this way the spatial coherence of spontaneous emission and amplified spontaneous emission is 
consistently described by the quantum Maxwell-Bloch equations. Typical features of the model have been illustrated 
by the spectra of gain and spontaneous emission. An example of the spatial coherence characteristics described 
by the quantum Maxwell-Bloch equations has been presented by analytically obtaining the spontaneous emission 
factor [3 and the far field distribution for the example of a broad area edge emitting laser. In general the quantum 
Maxwell-Bloch equations derived for edge emitting and vertical cavity surface emitting lasers provide a starting point 
for a detailed analysis of spatial coherence patterns in diverse semiconductor laser geometries such as broad area or 
ultra-low threshold lasers. 



APPENDIX A: STOCHASTIC SIMULATION OF MEASUREMENTS 

As explained at the end of section [v|, the field average of zero is an expression of our lack of knowledge about the 
actual physical field present in the laser device. This lack of knowledge could be removed by performing a measurement 
on the light field. For example, the phase information could be obtained by measuring the interference of the light field 
from the diode with light from a separate laser. A simple theoretical simulation of such a measurement is obtained by 
assuming that I(r; r') defines the variance of a Gaussian distribution of coherent states (Gaussian P representation). 
The measurement of the actual field E(r) may then be simulated by randomly selecting a coherent field from this 
Gaussian distribution. Since the field and the dipole density are correlated, the statistical selection of E(r) shifts the 
dipole density average p(r,k) from zero to 

_ /rfr'E(r)C(r';r,k) 

Because the field has been determined with a precision equal to the quantum limit, the correlations factorize and the 
new values of I(r; r') and C(r; r', k) are given by the respective products of E(r) and p(r, k). 

The drastic changes in the intensity distribution and in the dipole dynamics which may result from such a simulated 
measurement cause spatial holeburning and give rise to relaxation oscillations. By performing several stochastical 
simulations starting from the same initial probability distribution, the statistics of the carrier density fluctuations 
may be obtained. It is then possible to derive the correct linewidth enhancement factor a as well as the relaxation 
oscillation sidebands directly from the dynamics of the quantum Maxwell-Bloch equations. 
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APPENDIX B: LANGEVIN EQUATIONS 



An alternative approach to the problem of quantum noise is given by the Langevin equations ]25|,?]. Langevin 
equations simulate quantum noise by adding classical noise sources which reproduce the statistical properties given 
by the uncertainty relations. In particular, the exponential damping terms given by the cavity loss rate k and the 
dipole relaxation rate T must be compensated by a noise input maintaining the quantum fluctuations in the field and 
dipole densities. For quantum well lasers, these fluctuations are given by 



{£*(r)£ j (r)) = -5 ij 5(r-v') (Bla) 



( P ;(r„,k||) ft „(r'|, > k'|,)) = ^„*(r|, - r',,)^ 2 (f! q (fc,,, JV(r„)) • (*||,JV(r,|)) 

+ (1 - f eq (*|,,JV(r||))) ■ (1 - f h eq (k h N(v l{ )))) . (Bib) 

Note that the variables £i(r) and p(r y , k||) are classical quantities in this context. The Langevin equations presented 
in the following do not describe quantum mechanical coherent states. Instead, semiclassical light fields and dipole den- 
sities are calculated with a precision violating the uncertainty relations. While the statistics thus derived correspond 
to the quantum statistics, the theory applied is fully classical. 

For a VCSEL, the semiclassical Langevin equations corresponding to the quantum Maxwell-Bloch equations ( pOj 
read 

^iV(r) = £> am6 AJV(r) + j(r) - 7 7V(r) 



fVQ 



J d 2 k (j20*(r)pi(r, k) - £(r)# (r, k)) 



^ft(r s k) = -(r(k)+ifi(k))p i (r J k) 



r 



— £i(v) = -{m + i8uji)£i(v) + i— — A£j(r) 



-igo- 



I d 2 k ^( r ,k)+^^-^-^(r,k)j +qf(r,t) (B2c) 



4tt 2 . 

with the noise input terms qf (r, k, t) and qf (r, t) given by 

(qf *(r, k, t)<£(r', k, i + t)> = 5(r)r(k)* y 5(r - r')47r 2 a (/ e e g (fc, N(rj) ■ (fc, AT (r) 

+ (l - f e eq (k, N(r)) ) • (l - /* (k, N{v)) ) ) (B3a) 

(qf(r,t)qf(r>,t + T)) = 5(T) K S(r-r>). (B3b) 

These Langevin equations allow a stochastical simulation of quantum noise in a semiclassical context. However, it is 
necessary to perform a large number of runs before statistical results can be obtained. In particular, the unrealistic 
intensity fluctuations of the light field vacuum need to average out before the results are consistent with quantum 
theory. 
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FIG. 1. Schematic representation of the edge emitter geometry. The laser field is mostly confined to the plane of the quantum 
well and propagates along the y-axis. 

FIG. 2. Schematic representation of a typical VCSEL geometry. The laser filed propagates perpendicularly to the quantum 
well along the z-axis. The length of the optical resonator approximately corresponds to the wavelenth. 

FIG. 3. Spectra of gain and spontaneous emission for a quantum well structure of Q — 5 quantum wells given relative to 
G m ax = 2.5 x 10 7 cms _1 <5(r 2 — z ) and S max = 2.7 x 10 7 cm~ 2 . The five spectra shown are the gain at N — (go), the gain at 
N — 10 12 cm~ 2 (gi), the spontaneous emission at N = 10 12 cm~ 2 (si), the gain at N — 5 x 10 12 cm~ 2 (g5) and the spontaneous 
emission at N = 5 x 10 12 cm~ 2 (s5). 

FIG. 4. Carrier density dependence of the spontaneous emission factor f3 for three modes with frequencies above the band 
gap frequency given by (a) lo — 0, (b) ui = 0.5(r + n), and (c) uj — T + k. f3o = (3(u> = N = 0) is determined by the geometry 
of the laser. The carrier density is given in terms of the transition frequency at the Fermi edge flf. 

FIG. 5. Far field intensity distributions for R = 0.5, %ojq = 1.5 eV, 7i(F + k) = 8 meV and e r = 12. The density is pinned 
at Qf = 1.8805(r + k). Figure (a) shows the far field pattern for carrier densities of 0.05, 0.10 and 0.15 times pinning density, 
figure (b) shows the distribution for 0.25, 0.5. and 0.75 times pinning density and figure (c) shows the distribution for 0.90, 
0.95 and 0.99 times pinning density. The peaks appear at emission angles of ±15°. 
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